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RESUMEN 

Para el analisis no lineal de estructuras, es necesario conocer el diagrama 
momento- rotacion o el diagrama momento-curvatura, para las vigas y columnas 
de hormigon armado o de acero y para las diagonales rigidizadoras se necesita 
conocer el diagrama fuerza-desplazamiento. El ASCE 41 de 2011, proporciona 
tablas que facilitan el calculo de estos diagramas, las mismas que se presentan en 
este articulo, con ejemplos desarrollados, manualmente. 

Entre los ejemplos se presenta una estructura de 4 pisos de hormigon 
armado y una estructura de acero de 6 pisos; en las dos, se tiene disipadores de 
energia ADAS o TADAS, colocados sobre diagonales de acero en forma de “V” 
invertida; se presentan los diagramas momento curvatura en varios elementos de 
estas dos estructuras y manualmente se ilustra el calculo del diagrama fuerza- 
desplazamiento para las diagonales de acero, tanto para el caso de compresion 
como para el caso de traccion; de igual forma se obtiene manualmente el 
diagrama momento-curvatura para un disipador de energia. 

El calculo con las tablas del ASCE 41 implica evaluar algunos criterios y la 
mayor parte de las veces se debe interpolar con los valores de las tablas, de tal 
manera que es bastante laborioso hacerlo manualmente por lo que se presentan 
varios programas del sistema de computacion CEINCI-LAB que obtienen los 
diagramas indicados en el primer parrafo. 
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ABSTRACT 

For nonlinear analysis of structures, it is necessary to know the moment- 
rotation diagram or moment curvature diagram, for beams and columns of 
reinforced concrete or steel and for the stiffening diagonals, these needs to know 
the force-displacement diagram. The ASCE 41 of 2011 provides tables that 
facilitate the calculation of these diagrams, the same that those presented in this 
article, with examples developed manually. 

In the examples are presented a reinforced concrete structure of 4 floors 
and a steel structure of 6 floors, in the these structures, they have energy 
dissipators TADAS or ADAS mounted on steel bracing in the form of inverted "V"; 
the moment curvature diagrams are presented in various elements of these two 
structures and manually calculated the force-displacement diagram for steel 
diagonal, for both the case compression as in the case of traction, likewise 
manually, it is obtained the moment-curvature diagram for an energy dissipators. 

The calculation with the table of ASCE41 involves assessing certain criteria 
and most of the time must be interpolated with the values of the tables, so it is 
quite laborious to do it manually, so several programs are presented of computer 
system CEINCI-LAB. They get the diagrams indicated in the first paragraph. 



1. INTRODUCTION 

Para el analisis no lineal: estatico o dinamico, se necesita conocer las 
relaciones momento-curvatura o momento-rotacion, para el caso de flexion y flexo 
compresion; de estos diagramas se obtiene la rigidez de una seccion de acuerdo 
al nivel de cargas que gravita sobre la misma. (Aguiar et al., 2015, 1; Mora y 
Aguiar, 2015; Medina, 2012). 

En el sistema de computacion CEINCI-LAB se tienen los programas 
MCvIga y MCcolumna que determinan el diagrama momenta curvatura de vigas y 
columnas, respectivamente, utilizando un modelo de hormigon no confinado y un 
modelo elasto plasto para el acero. Ortiz y Freire (2013) 

Para hormigon confinado tambien se han desarrollado programas para 
encontrar el diagrama momenta curvatura, utilizando el modelo tri lineal para el 
comportamiento del acero. Aguiar (2002). En resumen el diagrama momenta 
curvatura depende de los modelos constitutivos del hormigon y del acero que se 
utilicen y para cualquiera de los casos con los que se trabaje se necesita contar 
con programas de ordenador. Borda (2007), Chioul et al. (2008), Rogac (2012), 
Lepe (2013). 

Ahora bien, el ASCE 41 presenta unas tablas que tienen un respaldo 
teorico y experimental, por medio de las cuales se puede encontrar en forma 
sencilla las relaciones momenta rotacion y a partir de estas inferir las relaciones 
momenta curvatura con las que se halla la rigidez de acuerdo en una seccion de 
acuerdo al estado de dano que tiene la misma. 





Diagrama momento-rotacion y momento-curvatura 
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Figura 1 Relaciones momento curvatura y modelo del ASCE 41 

A la izquierda de la figura 1, se presenta las relaciones momento 
curvatura, para un modelo que contempla tres zonas de trabajo del elemento; la 
elastica (sin dano) que es para cuando el momento actuante M < M Y \ siendo M Y 
el momento de fluencia; en esta zona la rigidez a flexion vale E I e . Donde E es el 
modulo de elasticidad del material; I e es el momento de inercia de la seccion 
gruesa, que se halla en base a la geometrta de la seccion; ast para una viga 
rectangular de seccion constante l e = b h 3 /12, donde b,h, son la base y altura de 
la seccion; de tal manera que la rigidez elastica se halla en base a la geometrta de 
la seccion transversal del elemento sin utilizar el ASCE 41 . De igual manera se 
halla el momento de fluencia M Y . 

La segunda zona es la de post fluencia, en esta se cumple M Y < M < 
M u , y la rigidez de la seccion es E / p ; para esta zona se requiere las propuestas 
del ASCE 41 pero antes de explicar su calculo y para terminar con las zonas se 
indica que la ultima es la residual que tiene un momento M R . 

El diagrama presentado a la izquierda de la figura 1 es para el caso de 
secciones asimetricas que tienen un comportamiento diferente si la parte inferior 
esta traccionada y la superior trabaja a compresion; o al reves. 

A la derecha de la figura 1 , se presenta la nomenclatura del ASCE 41 para 
hallar las relaciones momento rotacion. En el eje horizontal se indica 9 (giro) o 
desplazamiento A; y en el eje vertical Q/Q Y . Para el caso de flexion esta relacion 
sera M/M Y y lo interesante es que el valor para el punto B, es la unidad, de tal 
manera que el momento es el de fluencia M Y , que esta asociado a una rotacion e Y . 

A partir del punto M Y ,9 Y , se encuentra el punto U (grafico de la izquierda) 
que corresponde al punto C (grafico de la derecha) por medio de la variable a 
propuesta por ASCE 41. Se determina tambien el momento residual, M R , con la 
variable c; y la rotacion final en este segmento ~DE por medio de la variable b. Se 
destaca que en el hormigon el punto U esta asociado a una deformacion a la 
compresion que esta entre 0.003 y 0.004; en cambio para la zona residual esta 
deformacion se halla entre 0.004 y 0.008. 
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Las variables a,b,c, presentan el ASCE 41 de 2011 para elementos de 
acero y hormigon armado y se indican en este artlculo con ejemplos y programas 
de computacion del sistema CEINCI-LAB. 

Es una gran ventaja contar con las tablas del ASCE 41 para hallar las 
relaciones momento-rotacion, en forma relativamente sencilla pero se debe tener 
claro que en el modelo del ASCE 41 (figura 1) el punto C, se encuentra a la 
derecha del punto B, es decir la seccion tiene ductilidad por curvatura pero si la 
viga no tiene ductilidad debido a que la cuantfa de acero a traccion es similar o 
igual a la cuantfa balanceada, no va a existir el punto C. (Marin 1979). De tal 
manera que se recomienda un poco de precaucion, tratando de no estar en los 
Ifmites de los criterios que presentan las tablas. 

2. PUNTO DE FLUENCIA PARA VIGAS Y COLUMNAS DE 
HORMIGON 



Como se acaba de ver el diagrama momento rotacion se halla a partir del 
punto de fluencia Y, utilizando los coeficientes a,b,c, del ASCE 41. Para encontrar 
el punto Y se decidio trabajar con la propuesta de Y. Park (1985) que cuenta con 
un respaldo teorico y experimental del ensayo de 400 elementos de vigas y 
columnas de hormigon armado. Las ecuaciones de calculo se indican a 
continuacion. 



M Y = 0.5 fc b d 2 [(1 + P c ~ V) Vo + (2 - v) Pt + Ol ~ 2 p c )a c p' t ] 



d' U./b /£ c \ 

d ^ 1+ a v \£ n / 



0.75 /£ c x 0,7 



OCy Tj 0 



Pt = 



As fy 

b d 



, A 'sfy 

p ‘ = tdfJ 



bdf; 

8 C tpy d £y ^ £ u 



«C = ( 1 — Pc) -r- & <1 



«V = [l-05+ fe-1.05)^] 



0.03-1 (1 — k)d 



k = 



Kpt + p't y 



4 a 



\ + ( Pt+ Pc P't ) ~T ~(Pt+ Pt ) 



a,, 



1 

2 a y 



0.45 

C2_1 + (0.84 + p t ) 



Donde b,d,d', son la base, la altura efectiva y el recubrimiento de la 
armadura a compresion; / y , £ y , / c ', e 0 son: el Ifmite de fluencia del acero, la 
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deformacion de fluencia, la resistencia maxima del hormigon a compresion, y la 
deformacion asociada a la resistencia maxima del hormigon; P 0 es la carga axial 
que gravita sobre el elementos, para vigas es igual a cero; a c es la deformacion del 
acero a compresion que tiene que ser menor a la deformacion ultima e u que puede 
ser igual a 0.003 para diseno o 0.004 para analisis; A S ,A' S es el area de la 
armadura que trabaja a traccion y compresion respectivamente, para el caso de 
columnas se considera que A s = A’ s ; M Y ,0 Y , son el momento y curvatura en el 
punto de fluencia. 

Con el modelo de Park (1985) tambien se puede hallar el momento y 
curvatura para el punto ultimo (M U ,0 U ), en Aguiar (2003) se indica las ecuaciones 
de calculo para este punto. 



3. RELACIONES PARA VIGAS DE HORMIGON 



En las tres primeras columnas de la tabla 1, se indican los parametros que 
se deben evaluar, antes de encontrar los valores de a,b.c. que estan en las tres 
ultimas columnas. El primer parametro esta relacionado con las cuantfas de acero; 
el segundo define si el refuerzo transversal utilizado esta confinando o no a la viga 
y el tercero esta relacionado con el cortante actuante. A continuacion se indica el 
calculo de estos parametros. 

• Cuantfa a flexion 

_ A s f _ A' s _ 0.85 f' o ( 6120 \ 

P bd P bd Pb f y Pl \6l20 + f y ) 

/?! = 0.85 para f'c < 280 kg /cm 2 

Pt = 0.65 < 1.05 - /'c/1400 < 0.85 Para f'c > 280 kg /cm 2 

Las variables todavfa no definidas son: p,p',p b , que son la cuantfa del 
acero a traccion, a compresion y la cuantfa balanceada. Con estas 

cuantfas de determina la relacion: y en base a su valor se ingresa a 



• Refuerzo Transversal confinado, si cumple: 



s 




3 

V s>7 V 

4 



V s = 



A v * fy * d 

s 



Donde A v es el area del refuerzo transversal; f y es el Ifmite de fluencia; d 
es la altura efectiva; s es el espaciamiento de los estribos; es la 
capacidad al code de la armadura transversal; V es el cortante actuante. 






Cortante actuante 
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V 

K * d * 4JI 

Si se trabaja el cortante V en T.; las dimensiones de la viga o columna: 
b w ,d, en m., y la resistencia del hormigon en T/m 2 , se debe multiplicar el 
resultado por 1.1926 para pasar Pound/in. 



Tabla 1 Parametros y criterios de aceptacion para procedimientos no 



lineales en vigas de hormigon. 



i.i) Vigas controladas por flexion 


a 


b 


C 


P~P' 


Refuerzo 

Transversal 


V (pound) 








Pbal 


b w (in)d(iri)J f'c(psi ) 


< 0.0 


Confinado 


< 3 


0.025 


0.05 


0.2 


< 0.0 


Confinado 


> 6 


0.02 


0.04 


0.2 


> 0.5 


Confinado 


< 3 


0.02 


0.04 


0.2 


> 0.5 


Confinado 


> 6 


0.015 


0.02 


0.2 


< 0.0 


No Confinado 


< 3 


0.02 


0.03 


0.2 


< 0.0 


No Confinado 


> 6 


0.01 


0.015 


0.2 


> 0.5 


No Confinado 


< 3 


0.01 


0.015 


0.2 


> 0.5 


No Confinado 


> 6 


0.005 


0.01 


0.2 


i.ii) Vigas controladas por corte 








Espaciamiento de estribos < d/2 


0.0030 


0.02 


0.2 


Espaciamiento de estribos > d/2 


0.0030 


0.01 


0.2 



• Ejemplo 1 

Dada una viga cuadrada de 40 cm., de lado y doblemente armada como 
se indica en la figura 2; el recubrimiento es de 6 cm.; el refuerzo transversal esta 
constituido por varillas de 10 mm, espaciados cada 10 cm. El hormigon tiene un 
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f' c = 240 kg /cm 2 . y el acero es A 36. El modulo de elasticidad del hormigon es 
1500000 T/m , del acero 20430000 T/m 2 . Se pide: 



i) Determinar la rigidez axial EA, y rigidez a flexion El en el rango elastico. 

ii) Encontrar los parametros a, b, c, de la tabla 1 si el cortante V = 0 y el 
momento inicial M = 0, inicio de calculo. 

iii) Hallar el Momento y curvatura a nivel de fluencia. 



E01Omm @10 cm 



S? 




„ .40 

^ 



E c = 1'500.000 T/m 2 
f y = 25300 T/m 2 
E s = 20'430.000 T/m 2 
f u = 40800 T/m 2 



Figura 2 Geometrfa de una seccion de una viga y parametros del 



acero y concreto. 



• Solucion 

i) Rigidez axial y a flexion 

A = 0.40 * 40 = 0.16m 2 
0.40 * 40 3 

/ = — = 0.00213 m 4 

12 

La rigidez axial EA y a flexion El, valen: 

EA = 240000 T. El = 3200 T/m 



ii) Factores a, b, c, de tabla 1 . 

4s = 3016 mm = 6.03 As' = 3020 mm = 9.42 cm 2 



6.03 
40 * 34 



0.0038 



9.42 
40 * 34 



0.0059 



P~P' 



-0.4039 



p b = 0.0053 



Pb 
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A continuacion, se procede a determinar si la seccion se encuentra 
confinada, para lo cual la viga se encuentra armada con 1<t> 10 mm @10cm. 

A v — 2010 mm — 1.57 cm 2 

1.57 * 4200 * 34 

V s = — = 22.43 T 

5 10 

Al cumplir las dos condiciones presentadas anteriormente se puede 
determinar que la seccion se encuentra confinada. Por otra parte, al ser V = 0, 
dato del problema. Los valores a,b,c son los de la primera fila de la tabla 1 . 

a = 0.025 b = 0.050 c = 0.200 
iii) Momento y curvatura de fluencia 

Al aplicar las ecuaciones de Y. Park (1985) se halla que el momento y 
curvatura de fluencia valen: 



rad 

M y = 7.6679 Tm. 0 r = 0.0081 

m 

3.1 Longitud plastica 

Cuando se supera el momento de fluencia M Y se tiene una longitud 
plastica en el nudo inicial que se denomina L pa , y cuando se supera el momento 
de fluencia M Y en el nudo final se tiene la longitud L pb \ a la izquierda de la figura 3 
se presenta la nomenclatura indicada donde los momentos actuantes se han 
denominado en la grafica estos dos momentos son anti horarios. 





Figura 3 Diagrama de momentos y longitudes plasticas 

En la grafica de la derecha de la figura 3, el valor de M,- se lo ha colocado 
al lado del momento M u ya de que de esta forma se facilita el calculo de la longitud 
plastica. En efecto, de la relacion de triangulos semejantes, se tiene: 

Mi + Mj _ M Y + Mj 

L L Lp a 
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Al despejar L pa se encuentra: 

_ ( Mj — My \ 

pa yMj + M ; J 

Para el nudo final, se procede en forma similar y se halla: 

_ ( Mj ~ 

pb " ^ + Mj) 

Notese en las dos ecuaciones que el numerador debe ser siempre 
positivo, de tal manera que la longitud plastica sera positiva. 

3.2 Rotacion de fluencia 

Para calcular la rotacion de fluencia d Y es necesario recordar la forma 
como se obtiene la matriz de flexibilidad en una viga que se encuentra sobre un 
apoyo fijo en el nudo inicial y un rodillo en el nudo final. Aguiar (2014). 

Para encontrar la matriz de flexibilidad se aplica un momento unitario en el 
nudo inicial y se calculan los giros en el nudo inicial y final que son L/3EI, L/6EI, 
respectivamente, grafico superior izquierdo de figura 4; al aplicar el momento 
unitario en el nudo final se hallan los giros indicados en la grafica superior derecha 
de la figura 4. Se recuerda la convencion de signos que anti horario es positivo. 




Figura 4 Elementos de la matriz de flexibilidad de un elemento de seccion 

constante 



Ahora, cuando se aplican los dos momentos unitarios a la vez, al aplicar el 
principio de superposicion lineal se tienen los giros indicados en la parte inferior de 
la figura 4. 

Si en lugar de aplicar un momento unitario, se tiene un momento M Y en el 
nudo inicial, y un momento M Y en el nudo final; los giros que se encuentran son los 
de fluencia 9 Y ,6 Y . Nuevamente por el principio de superposicion, estos giros valen: 
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LM y L My 
3 El ~ 6 El 

L My LMy 

3 EI~ 6 El 




Figura 5 Rotaciones de fluencia en nudo inicial y nudo final 

La figura 5 considera que las dos secciones ingresan simultaneamente al 
rango no lineal (hipotesis conservadora) ya que la mayor parte de las veces solo 
una de las secciones ingresara al rango no lineal y la otra no. 

Para el caso particular en que M Y = M Y , la rotacion de fluencia es la 
recomendada por el ASCE 41, que resulta la siguiente para elementos de 
hormigon armado. 



e Y = 



LMy 

6 El 



Para elementos de acero el momento de fluencia M Y = zf y . Luego la 
ecuacion resulta. 



8y = 



Lzfy 

6 El 



Las variables todavla no indicadas, son: z es el momento estatico (z = 
/ y dA = Ac * brazo_s + At * brazoj)', f y es el esfuerzo de fluencia del acero. 



Las ecuaciones indicadas son para vigas. Ahora para columnas, el ASCE 
41 recomienda la siguiente ecuacion. 



8y = 




Donde P es la carga axial que gravita en la columna; P CL es la carga crltica 
de pandeo. De acuerdo al ASCE 41, la carga crftica se calcula de la siguiente 
manera. 



Se determina la relacion L/r. Donde L es la longitud del elemento; r es el 
radio de giro. Se debe verificar que: 
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L 

- < 108 
r 

Si la relacion L/r es menor a 108; el esfuerzo del material F cr es 12 (psi); 
caso contrario vale: 

1.4* 10 5 

F cr = r~2 — (.Psi) 

(f) 

Para pasar de (psi) a (T/m 2 ) se debe multiplicar por 703.069. Finalmente, 
la carga crltica de pandeo se halla con la siguiente ecuacion. 

P = AF cr 

Donde A es el area de la seccion transversal. 



3.3 Factor que mide efecto de corte 



El factor que mide el efecto de corte 0 para un elemento de seccion 
constante que trabaja en el rango elastico vale (Aguiar, 2014) 



0 = 



3 PEI 
GAL 2 



Las variables no definidas son: /? coeficiente de forma que vale 1.2 para 
secciones rectangu lares; G es el modulo de corte. Ahora bien, interesa calcular el 
factor 0 para un elemento de seccion constante, que ha ingresado al rango no 
lineal y tiene una longitud plastica L pa asociada a una rigidez El a , en el nudo 
inicial; en el nudo final la longitud plastica es L pb y la rigidez es El b , la longitud 
central que no ha ingresado al rango no lineal es: L — L pa — L pb , y tiene una 
rigidez EI 0 . Solo para el calculo del factor 0 se determinar una rigidez ponderada 
en funcion de las longitudes, denominada El m . 

pr _ Lpa Pi a + (L — L va — Lp EI q + Lp b EI b 
^ 'm ~ r 



Finalmente se determina un solo factor del efecto de corte, aplicable al 
nudo inicial, centra de luz y nudo final, denominado 0 m 



0m 



3P EJm 

GAL 2 



• Ejemplo 2 

La viga de 40/40 del Ejemplo 1, tiene una longitud de 6.0 m.; un momento 
de fluencia M Y = 7.6679 Tm., una curvatura de fluencia 0 F = 0.0081 — . En un 

m 
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instante de cargas dinamicas, el momento actuante en sus extremos es de 
Mi = Mj = 8 Tm., y el cortante actuante es V = 2.7 T. Se pide: 

i) Calcular la rotacion de fluencia d Y de acuerdo al ASCE 41 ; 

ii) La longitud plastica, considerando M Y = M Y . 

iii) El diagrama Momento rotacion con la formulacion ASCE 41 

iv) El diagrama Momento curvatura de acuerdo al ASCE 41 



• Solucion 

i) Rotacion de fluencia 



_ 6 * 7.6679 
6y ~—* 3200 



0.0024 



ii) Longitud plastica 

Los momentos actuantes en los extremos son anti horarios y valen 8 Tm., 
(dato del ejemplo), los momentos de fluencia valen 7.6679 Tm, como se 
indica en la figura 6. 




Figura 6 Momentos actuantes y de fluencia para calcular longitud plastica 



/8.00 - 7.6679\ 

= — — — *6 = 0.1245 m 

V 8.00 + 8.00 / 



iii) Relacion momento rotacion 

En el Ejemplo 1, se encontro que ( p - p')/p b es menor a cero; que la viga 
esta confinada; y ahora, en base al cortante actuante de V = 2.7 T. Se 
demuestra que la relacion V/(b d JfJj es menor a 3. 

V 2.7 Pound 

— = == * 1.1926 = 1.5283 

b d yffj 0.4 * 0.34 * V240 in 



Luego de la tabla 1 se encuentra: 
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a = 0.025 b = 0.050 c = 0.200 

Con estos valores se halla el diagrama momento rotacion de acuerdo al 
ASCE 41 

6 V = 6 Y + a = 0.0024 + 0.025 = 0.0274 rad 
Mu=M y + 0.05 El (i 9a ~ 0 Y ) 

Mjj = 7.6679 + 0.05 * 3200 (0.0274 - 0.0024) = 11.6679 T m 
M R — c M y — 0.20 * 7.6679 = 1.5336 T m 
0 R = 8 Y + b = 0.0024 + 0.05 = 0.0524 rad 
El diagrama momento rotacion se indica en la figura 7. 




iv) Diagrama momento curvatura 



rad 



0 Y = 0.0081 



m 
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dr, 0.0274 rad 

0t/ = + -T = 0.0081 + — — = 0.2282 

L p 0.1245 m 

6 r 0.0524 rad 

0B = 0y+ _ = 0.008! + — = 0.4290 — 



En la figura 8 se indica el diagrama momento curvatura. 




4. RELACIONES PARA COLUMNAS DE HORMIGON 

En forma similar a las vigas, hay tres controles que se deben realizar, antes de 
ingresar a la tabla 2 para encontrar los coeficientes a,b,c con los cuales se halla el 
diagrama momento rotacion. El primer control es calcular la relacion: 

P 

Ag fc 

Donde P es la carga axial; A g es el area de la seccion transversal de la 
columna y / c ' es la resistencia maxima del hormigon a compresion. El segundo 
control tiene que ver con el confinamiento y el tercero con el cortante actuante; 
estos dos controles se realizan en forma similar a la indicada para vigas. 

Si la columna esta controlada por code, ultima fila de la tabla 2, no se 
puede obtener los coeficientes a,b,c. 

Tabla 2 Parametros y criterios de aceptacion para procedimientos no lineales 



en columnas de hormigon armado. 



ii.i) Columnas controladas por flexion 


a 


b 


c 



